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We present the results of the lattice QCD calculation of the average up-down and strange quark 
masses and of the light meson pseudoscalar decay constants, recently performed with Nf = 2 
dynamical fermions by the ETM Collaboration. The simulation is carried out at a single value 
of the lattice spacing with the twisted mass fermionic action at maximal twist, which guaran- 
tees automatic 6 (a) -improvement of the physical quantities. Quark masses are renormalized 
by implementing the non perturbative RI-MOM renormalization procedure. Our results for the 
light quark masses are mjjf (2 GeV) = 3.85 ± 0.12 ± 0.40 MeV, mp(2 GeV) = 105 ± 3 ± 9 
MeV and m s /m ud = 27.3 ±0.3 ± 1.2. We also obtain f K = 161.7 ± 1.2 ± 3.1 MeV and the ra- 
tio fic/fa = 1 -227 ± 0.009 ± 0.024. From this ratio, by using the experimental determination of 
T(K — » pLV^(y)) /Y(n — > M^V(y)) an d tne average value of \V uc i \ from nuclear beta decays, we ob- 
tain |V MJ | = 0.2192(5) (45), in agreement with the determination from Kb decays and the unitarity 
constraint. 
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Figure 1: Effective masses of pseudoscalar mesons with jis = jii = jJ-2, as a function of the time. 

1. Introduction 

We present our recent determination ^ of the light quark masses (strange quark mass m s and 
average up-down quark mass of the kaon pseudoscalar decay constant fx, and of the ratio 
Ik/ fn- m order to investigate the properties of the K meson, we have simulated the theory with 
Nf = 2 degenerate dynamical and two valence quarks, considering a partially quenched setup with 
the valence quark masses /ii and /i 2 different between each other and from the sea quark mass ii s . 

The calculation is based on a set of (ETMC) gauge field configurations generated with the tree- 
level improved Symanzik gauge action at j8 = 3.9, corresponding to a = 0.087(1) fm (a -1 ~ 2.3 
GeV) |0], and the twisted mass fermionic action at maximal twist. We have simulated 5 val- 
ues of the bare sea quark mass, ati s = {0.0040, 0.0064, 0.0085, 0.0100, 0.0150} and 8 values, 
aiiy 2 = {0.0040, 0.0064, 0.0085, 0.0100, 0.0150, 0.0220, 0.0270, 0.0320}, for the valence quark 
mass. The first five masses, equal to the sea masses, lie in the range 1 /6m s ^ jUi,2 ^ 2/3m s , being 
m s the physical strange quark mass, while the heaviest three are around the strange quark mass. 

At each value of the sea quark mass we have computed the two-point correlation functions of 
charged pseudoscalar mesons, on a set of 240 independent gauge field configurations, separated by 
20 HMC trajectories one from the other. To improve the statistical accuracy, we have evaluated the 
meson correlators using a stochastic method with a Z (2) -noise to include all spatial sources |]]. 
Statistical errors on meson masses and decay constants are evaluated using the jackknife procedure, 
while those on the fit results, based on data obtained at different sea quark masses, are evaluated 
using a bootstrap procedure. Further details on the numerical simulation can be found in refs. [|I|, [5]]. 

The use of twisted mass fermions presents several advantages ^ : i) the pseudoscalar meson 
masses and decay constants are automatically improved at G{a); ii) at maximal twist, the physical 
quark mass is directly related to the twisted mass parameter of the action, and it is subject only to 
multiplicative renormalization; iii) the determination of the pseudoscalar decay constant does not 
require the introduction of any renormalization constant, and it is based on the relation 

M ps 

The meson mass Mp$ and the matrix element |(0|P 1 (0)|P)| have been extracted from a fit of the 
two-point pseudoscalar correlation function in the time interval t/a G [10,21]. In order to illustrate 
the quality of the data, we show in fig. [j] the effective masses of pseudoscalar mesons, as a function 
of the time, in the degenerate cases fa = Hi = 
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2. Quark mass dependence of pseudoscalar meson masses and decay constants 

The determination of the physical properties of K mesons requires to study the corresponding 
observables over a large range of masses, from the physical strange quark down to the light up- 
down quark. In ref. we have studied the quark mass dependence of pseudoscalar meson masses 
and decay constants by considering two different functional forms: i) the dependence predicted by 
continuum partially quenched chiral perturbation theory (PQChPT), ii) a polynomial dependence. 

PQChPT fits: Within PQChPT we have considered the full next-to-leading order (NLO) ex- 
pressions with the addition of the local NNLO contributions, i.e. terms quadratic in the quark 
masses, which turn out to be needed for a good description of the data up to the region of the 
strange quark. The PQChPT predictions ^ can be written as 



fr(&-fr)ln2fr &(&-&)ln2fr 



+a v £i 2 + a s £s + av v £12 + ass £j + a vs £12 £, s + a V D %m 2 ) > ( 2 - 1) 

fps(Ms,Hi,pi) = /■ 

+(b v + 1/2) £ 12 + (b s - 1/2) & + v ^i 2 2 + 6ss ^ 2 + §12 & + ^VD §m 2 ] , 

where & = 2B Hi/(4Kf) 2 , £ 7 = B (^ + Hj)/(4nf) 2 and % Dij = B (jm- fij) / (4iuf) 2 . The param- 
eters Bo and / are the LO low energy constants (LECs) , whereas ay, as, by and bs are related to 
the NLO LECs by ay = 4a% — 2a^, as = &GC6 — 4ot4, by = CC5, bs = 2a%. The quadratic mass terms 



1 - Inl^s - § 2 5 ln2& s + _^ } ' 2 In (| ) + 



in eq. (2. 1 ) represent the local NNLO contributions. The chiral logarithms, also known at two loops 
in the partially quenched theory ||], involve a larger number of NLO LECs whose values cannot 
be fixed from phenomenology in the Nf = 2 theory. Introducing their contribution would increase 
significantly the number of free parameters, thus limiting the predictive power of the calculation. 

Aiming at a percent precision, the impact of finite size corrections cannot be neglected in our 
study, where the lattice spatial extension is L = 24a ~ 2.2 fm and MpsL > 3.2. Since we have not 
performed yet a systematic study on different lattice volumes, we have estimated the finite size 
effects by including in the fits the corrections predicted by one- loop PQChPT [g] (for their explicit 
expressions see ref. [[I]]). 

Polynomial fits: The inclusion of the local NNLO contributions in the PQChPT predictions 



of eq. (2. 1 ) is required by the observation that the pure NLO predictions are not accurate enough to 
describe the quark mass dependence of pseudoscalar meson masses and decay constants up to the 
strange quark region. Not having considered the full NNLO chiral predictions, we have evaluated 
the associated systematic uncertainty, considering as an alternative description a simple polynomial 
dependence on the quark masses, for both pseudoscalar meson masses and decay constants: 

M^ s {jis,nun 2 ) =B (n 1 + fi 2 )- [l+avZ n + a s £ s + a vv tf 2 + assZs+ a vsZnZs + avD£m2] , 
fps(ps,Huli2) =/■ [1 + (by + 1/2)62 + (b s -l/2)&+ byytf 2 + bssZf + bvsZnZs + bvDZm2] ■ 

(2.2) 

The differences between the results obtained by performing either chiral or polynomial fits have 
been included in the final estimates of the systematic errors. 

'The pseudoscalar decay constant / is normalised such that f % = 130.7 MeV at the physical pion mass. 
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3. Chiral extrapolations 

The input data of our analysis [P are the lattice results for the pseudoscalar meson masses and 
decay constants obtained at each value of the sea quark mass, with both degenerate and non degen- 
erate valence quarks. We have excluded from the fits the heaviest mesons having both the valence 
quark masses in the strange mass region, i.e. with a/J.i2 = {0.0220, 0.0270, 0.0320}, considering 
therefore 150 combinations of quark masses. The number of free parameters in the combined fit of 
Mp S and fps is 14, but a first analysis shows that some of them (from 1 to 5 depending on the fit) 
are compatible with zero within one standard deviation, and are kept fixed to zero. 

In order to extrapolate the pseudoscalar meson masses and decay constants to the points cor- 
responding to the physical pion and kaon, we have considered three different fits: 

• Polynomial fit: a polynomial dependence on the quark masses is assumed for pseudoscalar 



meson masses and decay constants, according to eq. ( |2.2[ ). 



PQChPT fit: pseudoscalar meson masses and decay constants are fitted according to the 



PQChPT predictions of eq. (2.1 ) including the finite volume corrections derived in ref. 



Constrained PQChPT fit: this fit, denoted as C-PQChPT in the following, deserves a more 



detailed explanation. The main uncertainty in using eqs. (gjj) and ( |2.2| ) to describe the 
quark mass dependence of Mp S and fp$ is related to the extrapolation toward the physical 
up-down quark mass. On the other hand, we have shown in ref. |Q] that pure NLO ChPT, 
with the inclusion of finite volume corrections, is sufficiently accurate in describing the lattice 
pseudoscalar meson masses and decay constants when the analysis is restricted to our lightest 
four quark masses in the unitary setup (i.e. /ii = \i% = /I5). In order to take advantage of this 
information, when performing the C-PQChPT fit we first determine the LO parameters So 
and / and the NLO combinations ay + as and by + b$ from a fit based on pure NLO ChPT 
performed on the lightest four unitary points. By using these constraints, the other parameters 
entering the chiral expansions of Mp S and fp$ are then obtained from a fit to eq. (2A) over 



the non unitary points. For consistency with the previous unitary fit, we exclude also in this 
case from the analysis the data at the highest value of sea quark mass, a/15 = 0.0150. 
We find that, though the quality of the fit is better in the polynomial case, all three analyses provide 
a good description of the lattice data, in the whole region of masses explored in the simulation, 
once the terms quadratic in the quark masses are taken into account. 

A potential problem in the partially quenched theory is the divergence of the chiral logarithms 



in the limit in which the light valence quark mass goes to zero at fixed sea quark mass (see eq. (2.1 )). 
This divergence does not affect the extrapolation of the lattice results to the physical point, since 
sea and the light valence quark masses are degenerate in this case. However, in order to verify that 
this unphysical behaviour of the partially quenched chiral logarithms does not modify the result 
of the extrapolation, we have repeated the analysis restricting both the polynomial and the chiral 
fits to the 30 quark mass combinations (26 in the case of the C-PQChPT fit) that, satisfying the 
constraint /I2 > /li = pis, are not affected by dangerous chiral logarithms. The comparison between 
the results obtained by considering the two different sets of quark masses is reassuring, as it shows 
that the effects of potentially divergent chiral logarithms are well under control in our analysis. 

The mass dependence of the pseudoscalar meson masses and decay constants is illustrated 
in fig. ^ where lattice data are compared with the results of the polynomial, PQChPT and C- 
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Figure 2: Lattice results for a 2 Mp S / j(a/Xi +0/^2) (left) and a fp$ (right) as a function of the valence quark 
mass with ajj,\ = apis- The solid, dashed and dotted curves represent the results of the three fits. 



PQChPT fits. We have shown in the plots the cases in which one of the valence quark mass (/ii) 
is equal to the sea quark mass, and the results are presented as a function of the second valence 
quark mass (112)- The points corresponding to the physical pion and kaon are thus obtained by 
extrapolating/interpolating the results shown in fig. || to the limits ii\ — > m uc i and ii 2 — > m s . 

To investigate the impact of finite volume corrections we have compared, for the pseudoscalar 
meson masses and decay constants, the PQChPT fits obtained with or without including these 
corrections. The differences turn out to be small |jl|]; however, to better quantify the systematic 
error due to finite size effects, we plan to extend our analysis on lattices with different spatial sizes. 



By having determined the fit parameters, we have then extrapolated eqs. (2.1) and (2.2) to the 
physical pion and kaon, as follows. We have first used the experimental values of the ratios M K /f K 
and Mk/M k to determine the average up-down and the strange quark mass respectively. Once these 
masses have been determined, we have used again eqs. ( |2T| ) and (12) to compute the values of the 
pion and kaon decay constants as well as their ratio fx/ fit- 



4. Physical results 

In order to convert into physical units the results obtained for the strange quark mass and 
the kaon decay constants we have fixed the scale within each analysis (polynomial, PQChPT and 
C-PQChPT fits) by using f n as physical input. In the case of the polynomial and PQChPT fits 
we conservatively introduce for the dimensionful quantities a 6% and 3% of systematic error to 
take into account the different scale estimate derived in the analysis over the lightest four unitary 
points %%. 

The determination of the physical strange and up-down quark masses also requires implement- 
ing a renormalization procedure. The relation between the bare twisted mass at maximal twist, \i q , 
and the renormalized quark mass, m q , is given by m q {iiR) = Z m {g 2 ,oiiir) pL q {a), where iir is the 
renormalization scale, conventionally fixed to 2 GeV for the light quarks. Z m is the inverse of 
the flavour non-singlet pseudoscalar density renormalization constant, Z m = Zp l . We have used 
the non-perturbative RI-MOM determination of Z P , which gives ZjF~ MOM (l/<?) = 0.39(1) (2) at 



j8 = 3.9 [10], and converted the result to the MS scheme at the scale iir = 2 GeV by using renor- 



malization group improved continuum perturbation theory at the N LO [11]. 
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Fit 


m^F (MeV) 


mf 5 (MeV) 


m s /m ud 


h (MeV) 


fn 1 fit 


Polynomial 
PQChPT 
C-PQChPT 


4.07(9)(33) 
3.82(15)(25) 
3.74(13)(21) 


109(2)(9) 
107(3)(7) 
102(3)(6) 


26.7(2)(0) 
27.9(2)(0) 
27.4(3X0) 


158.7(11)(89) 
160.2(15)(54) 
161.8(10)(0) 


1.214(8X0) 
1.225(H)(0) 
1.238(7)(0) 



Table 1: Results for the light quark masses and pseudoscalar decay constants, in physical units, from 
the polynomial, PQChPT and C-PQChPT fits, analysing only the combinations of quark masses satisfy- 
ing fa > Ml = Ms- The quoted errors are statistical (first) and systematic (second), the latter coming from 
the uncertainties in the determination of the lattice scale and of the quark mass renormalization constant. 



In table |] we collect the results for the light quark masses and pseudoscalar decay constants, in 
physical units, and for the ratios m s /m uc i and fx/ fit, as obtained from the polynomial, PQChPT and 
C-PQChPT fits. To be conservative, we consider the results obtained from the analysis of the quark 
mass combinations satisfying the constraint jj, 2 > H\ = jls which, though being affected by larger 
statistical errors, are safe from the effects of the potentially divergent chiral logarithms. In table |] 
we quote as a systematic error within each fit the uncertainty associated with the determination of 
the lattice spacing and of the quark mass renormalization constant. 

In order to derive our final estimates for the quark masses and decay constants, we perform 
a weighted average of the results of the three analyses presented in table [I] and conservatively 
include the whole spread among them in the systematic uncertainty. In this way, we obtain as our 
final estimates of the light quark masses the results 

m]^(2 GeV) = 3.85 ±0.12 ±0.40 MeV , m^(2 GeV) = 105 ± 3 ±9 MeV, (4.1) 

and the ratio 

ms/mud = 27.3 ±0.3 ±1.2, (4.2) 

where the first error is statistical and the second systematic. For the kaon decay constant and the 
ratio fx/ fit we obtain the accurate determinations 

fx = 161.7 ±1.2 ±3.1 MeV , f K /f x = 1.227 ±0.009 ± 0.024. (4.3) 

An interesting comparison of our results for the strange quark mass and the ratio fx/ fit with 
other lattice QCD determinations is illustrated in fig. [3| (see ref. ^ for the full list of references). 

An important finding of our analysis [jl|] is that the use of non-perturbative renormaliza- 
tion turns out to play a crucial role in the determination of the quark masses. The estimate 
ZjP _M0M (l/a) = 0.39(1)(2) obtained with the RI-MOM method is in fact significantly smaller 
than the prediction Zf PT (l/a) ~ 0.57(5) given by one-loop boosted perturbation theory (in the 
same RI-MOM renormalization scheme) [[K]]. Had we used the perturbative estimate of Zp we 
would have obtained m^(2 GeV) = 2.63 ±0.08 ±0.36 MeV and mp(2 GeV) = 72 ±2 ±9 MeV. 
As shown in fig. [| (left), our prediction for the strange quark mass in eq. < ]4- 1[ > is in good agree- 
ment with other determinations based on a non-perturbative evaluation of the mass renormalization 
constant. The non-perturbative renormalization method, therefore, is found to have an important 
impact that can be even larger than the quenching effect and that should be kept in mind, particu- 
larly when combining the lattice results to produce the quark mass final averages. 
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Figure 3: Lattice QCD determinations of the strange quark mass (left) and of the ratio /k/ fn (right) 
obtained from simulations with Nf = 2 and Nf = 2 + 1 dynamical fermions. The results are also compared 
with the PDG 06 averages [12] and, for fx/ fx, with the average from theK^ determination ofV us /|7T 



Our result for the ratio fx/fn can be combined with the experimental measurement of F(K — > 
pLV tJL (y))/T(% — > HVfi(y)) [0] to get a determination of the ratio |K.v|/|Kd| [f^l- We obtain 
iKil/lK^I = 0.2251 (5) (47), where the first error is the experimental one and the second is the 
theory error coming from the uncertainty on fn/ fx- It yields, combined with the determination 
\V ud \ = 0.97377(27) [0] from nuclear beta decays, the estimate \V US \ = 0.2192(5)(45), in agree- 
ment with the value extracted from decays, \V US \ = 0.2255(19) [[l3|], and leads to the constraint 
due to the unitarity oftheCKM matrix \V ud \ 2 + \V US \ 2 + \V uh \ 2 - I = (-3.7 ± 2.0) • 10" 3 . 
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